Unclassified 

SECURITY  CLASSIFICA 


1a  REPORT  SECURITY  »  » 

Unclassified 


2a.  SECURITY  CLASSIFICATION  AUTHORITY 


AD-A160  225 


TATION  PAGE 

>.  RESTRICTIVE  MARKINGS 


3.  OlSTR I SUTlON/AVAl LABILITY  OF  REPORT 


Approved  for  public  release;  distribution 
unlimited 


A.  PERFORMING  ORGANIZATION  REPORT  NUMSERIS) 


S.  MONITORING  ORGANIZATION  REPORT  NUMBERIS) 

AFOSR- TR-  8  5-0  781 


6a  NAME  OF  PERFORMING  ORGANIZATION  Sb.  OFFICE  SYMBOL  7a  NAME  OF  MONITORING  ORGANIZATION 
Division  of  Engineering  at  applicable) 

Brown  University  AFOSR 


•a  ADDRESS  (City.  Slat*  and  ZIP  Cod* I 

Providence,  RI  02912 


7b.  ADDRESS  (City.  Slam  mad  ZIP  Coda) 

Bldg.  410 

Bolling  AFB,  D.C.  20332-6448 


.  NAME  OF  FUNOING/SPONSORING 
ORGANIZATION 

AFOSR 


.  ADDRESS  (City.  Simla  aid  ZIP  Coda) 

Bldg.  410 

Bolling  AFB,  D.C.  20332-6448 


11.  TITLE  (Include  Saeurit y  damnification) 

On  Structure  Determination  for  Polynomial  It 


12.  PERSONAL  AUTHORIS) 

A.  E.  Pearson 


13a  type  op  report 


Reprint 


S.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 

AFOSR-82-0230 


10.  SOURCE  OF  FUNDING  NOS. 


WORK  UNIT 
NO. 


1«.  DATE  OF  REPORT  (Yr  . 

13  March  1985 


IB.  page  COUNT 


COSAT I  COOES 


FIELD 


XXXX 


2323322222 


IB  SUBJECT  TERMS  (Continue  on  reverie  if  aaeammry  and  IdanUfy  by  block  number) 

Structure  determination,  polynomial  input-output, 
Fourier  series  coefficients,  FFT  algorithm 


If.  ABSTRACT  (Continue  on  reverie  if  neceaamry  and  idantify  by  block  number) 

>The  problem  of  structure  determination  for  a  deterministic  class  of  polynomial  input- 
output  differential  systems  is  formulated  as  a  minimum  norm-discrete  time  optimal  control 
problem.  The  order  of  the  differential  equation  and  the  degrees  of  the  polynomials 
involving  the  input-output  variables  play  the  role  of  multiple  discrete-timees  while  the 
coefficient  parameters  play  the  role  of  a  discrete  control  variable.  The  bas^.s  of  the 
parameter  identification  techniques  is  Shinbrot's  method  of  moment  functionals  using 
linear  combinations  of  commensurable  sinusoids  as  the  modulating  functions.  Given  the 
system  input-output  data  on  a  finite  time  interval,  the  underlying  computations  involve 
calculating  a  finite  set  of  Fourier  series  coefficients  of  moments  formed  from  the 
data,  which  can  be  efficiently  carried  out  via  an  FFT  algorithm,  followed  by  a 
sequence  of  singularity  tests  performed  on  a  controllability  type  Cram  determinant  that,^ 
arises  for  the  formulation.  0  g 

C~~- — 


20.  dibtribution/availability  of  abstract 

UNCLASSIFIED/UNLIMITEO  G3  SAME  AS  RPT.  B  OTIC  USERS  □ 


22a  NAME  OF  RESPONSIBLE  INDIVIDUAL 


Dr.  Marc  Q.  Jacobs 


DO  FORM  1473, 83  APR 


22b.  TELEPHONE  NUMBER 
(Inatude  A  too  Coda) 

(202)767-4940 


EDITION  OF  1  JAN  73  IS  OBSOLETE. 


Unclassified _ 

SECURITY  CLASSIFICATION  of  THIS  PAGE 


>  JS 


afosrtr 


0  781 


O  U 


ON  STRUCTURE  DETERMINATION  FOR 
POLYNOMIAL  INPUT-OUTPUT  DIFFERENTIAL  SYSTEMS1 


Accession  For 

NTIS  GRAJfcl 
DTIC  TAB 
Unannounced 
Justification- 


A.  E.  Pearson 
Division  of  Engineering 
Brown  University 
Providence,  RI  02912 


By - - - 

Distribution/ 
Availability  Codes 
lAvall  and/or  ~ 
Jist  I  Special 


ABSTRACT 


The  probiem  of  structure  determination  for  a  deterministic  class  of  polynomial  input-output 
differential  systems  is  formulated  as  a  minimum  norm-discrete  time  optimal  control  problem. 
The  order  of  the  differential  equation  and  the  degrees  of  the  polynomials  involving  the  input- 
output  variables  play  the  role  of  multiple  discrete-times  while  the  coefficient  parameters  play 
the  role  of  a  discrete  control  variable.  The  basis  of  the  parameter  identification  technique  is 
Shinbrot’s  method  of  moment  functionals  using  linear  combinations  of  commensurable  sinusoids 
as  the  modulating  functions.  Given  the  system  input-output  data  on  a  finite  time  interval,  the 
underlying  computations  involve  calculating  a  finite  set  of  Fourier  series  coefficients  of  moments 
formed  from  the  data,  which  can  be  efficiently  carried  out  via  an  FFT  algorithm,  followed  by  a 
sequence  of  singularity  tests  performed  on  a  controllability  type  Gram  determinant  that  arises 
from  the  formulation. 
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1.  INTRODUCTION 

The  problem  of  structure  determination  has  been  considered  for  the  class  of  linear  systems, 
both  stochastic  and  deterministic,  in  determining  such  structural  parameters  as  the 
controllability /observability  indexes  [1-3].  However,  within  the  context  of  system  identification, 
relatively  little  attention  appears  to  have  been  given  to  any  class  of  nonlinear  systems.  One 
such  class  is  considered  here,  namely  those  single  input-single  output  deterministic  systems  that 
can  be  described  by  a  polynomial  input-output  differential  equation  of  finite  order  involving 
finite  degree  polynomials  of  the  input-output  variables.  Examples  of  this  class  would  include 
linear  feedback  control  systems  with  zero  memory  polynomial  type  nonlinearities  inserted  in  the 
forward  and  feedback  loops,  and  certain  classical  nonlinear  differential  equation  models  such  as 
the  Van  der  Pol  and  Duffing  equations.  The  major  assumptions  are  that  the  input-output  data 
is  noise-free,  or  can  be  made  essentially  noise-free  by  appropriate  filtering,  and  that  upper 
bounds  are  known  for  the  degrees  of  the  polynomials  and  the  order  of  the  differentia]  equation. 
Although  the  assumption  of  noise-free  data  is  idealistic,  the  technique  utilizes  Shinbrot’s  method 
of  moment  functionals  [4]  with  linear  combinations  of  commensurable  sinusoids  as  the  modulat¬ 
ing  function  basis,  thus  facilitating  a  certain  measure  of  selective  filtering  on  the  data  as  a  result 
of  the  need  to  calculate  only  a  finite  set  of  Fourier  series  coefficients  while  avoiding  the  necessity 
to  estimate  unknown  initial  or  boundary  conditions.  Previous  developments  of  this  approach 
have  focused  on  linear  differential  systems  [5,6],  a  class  of  bilinear  systems  [7],  and  polynomial 
input-output  differential  systems  of  the  type  to  be  considered  here  [8],  all  of  presumed  known 
order  and  structure  as  far  as  the  differential  equation  model  is  concerned. 

Following  the  formulation  of  the  problem  as  a  triple  indexed  discrete  system,  a  solution  is 
proposed  using  minimum  norm-time  optimal  control  ideas,  and  some  discussion  is  included 
regarding  noise  and  bandwidth  considerations. 

2.  FORMULATION 

Let  [« (f  ),y  (f )]  denote  the  input-output  pair  observed  free  of  significant  measurement 
noise  on  a  finite  time  interval  [0,T],  and  let  p  denote  the  differential  operator  d/dt  so  that 
p2=rf2/ rft2,  etc.  The  class  of  polynomial  differential  equations  relating  u(f )  and  y(t)  is  taken 
of  the  form 

SSL  [»(<)]*  -  o  (i) 

I  =0/  =0*  —0 

0<»<7\  «j(0,0)=0,  » =0,l,2..n 

where  the  a,(j,k)  represent  parameters  which  are  to  be  determined  along  with  the  structural 
parameters  (n,m,/),  the  latter  consisting  of  the  order  of  the  differential  equation  n  and  the 
highest  degrees  (ro,f)  of  the  polynomials  involving  the  input-output  variables.2  The  choice  of 
the  time  interval  length  T  will  be  discussed  in  Section  4. 

The  actual  system  order  n0  and  the  degrees  of  the  polynomials  (m0,f0)  are  assumed  to 
satisfy  the  set  membership  expressions 

n06{l,2..n},  m0e{l,2.  m},  /„  €  {1,2.7}  (2) 

where  (n  ,m  ,T)  represent  the  previously  mentioned  a  priori  upper  bounds  on  the  structural 

2  Existence  and  uniqueness  of  solutions  to  (1)  on  [0,  T\  is  tacitly  assumed  for  whatever  values  are  as¬ 
signed  to  the  «,( y  ,k  ),  given  that  the  data  is  bounded  and  piecewise  continuous. 


parameters  ( n  ,m  ,1 ).  In  order  to  focus  on  determining  (n0,m0,/0),  the  model  (1)  will  be  re¬ 
parametrized  as  follows: 


n  I 


0  -»(«)+  E  E«i(*)i>i",+<" [»(<)]*  +  E  E/>i(i)p''X<)]’  (3) 

•  =14=1  1=1/ =1 

+  E  E  E-K(y.*)i>i-,[«(<)l'l»(01* 

i  —Ij  =14  =1 

where  ,  represents  the  Kronecker  delta.3  The  presence  of  S * ,  serves  the  dual  purpose  of  assur¬ 
ing  that  the  highest  derivative  term  appears  linearly  as  pn  y{t )  and  that  causality  is  preserved 
in  the  input-output  relation.  In  turn,  this  assures  that  (3)  has  a  normal  form  state  space  reali¬ 
zation  of  the  form 

*(0  —  8(*(0.«(0) 

f(0  -  c'*(0 

where  e  is  a  column  vector,  prime  denotes  transpose,  and  as  pointed  out  in  [8]  the  nonlinear 
terms  in  g(*,u)  are  polynomials  in  the  pair  (c'z,u).  Normalizing  the  coefficient  a0(0,l)  to 
unity  in  passing  from  (1)  to  (3)  is  somewhat  arbitrary  but  can  be  seen  to  be  valid  if  the  transfer 
function  for  the  linearized  differential  equation,  i.e.,  linearized  about  u  =y  =0,  does  not  contain 
a  pole  at  the  origin. 

Shinbrot’s  method  of  moment  functionals  facilitates  the  conversion  of  a  linear  differential 
relation,  such  as  (1)  or  (3),  to  an  algebraic  equation  in  the  coefficient  parameters  via  the  use  of 
so-called  modulating  functions.  The  resulting  algebraic  equation  is  characterized  by  functionals 
on  the  data.  A  motivating  factor  for  the  developments  in  [5-8]  has  been  the  specification  of 
modulating  functions  comprised  of  linear  combinations  of  the  commensurable  sinusoids 
{sinAk;0t  ,cos ku>0t  },  k  —0,1  ..L ,  u0=2w/T,  because  the  functionals  on  the  data  that  result  from 
the  conversion  are  linearly  related  to  the  first  L  Fourier  series  coefficients  of  moments  formed 
from  the  data  which  can  be  calculated  to  any  desired  degree  of  accuracy  by  a  sufficiently  high 
order  discrete  Fourier  transform  (DFT).  In  turn,  these  DFT's  can  be  efficiently  evaluated  by  a 
Fast  Fourier  Transform  (FFT)  algorithm.  In  order  to  apply  these  relations  to  the  problem  at 
hand,  let  f(f )  denote  the  (2 L  +1)  dimensional  column  vector  of  sinusoids: 


f(t)  =  Col  ^l;cosw0f  ,sinw0*  ;cos2u/0f  ,sin2w0f;  •  .  ;cosL  ,sin£ur0tj 


(4) 


0<f<r,  w0  =  2  jt/  t. 


Following  the  procedure  outlined  in  [5,6]  and  given  the  upper  bound  integer  a  for  the  system 
differential  equation  order,  let  C  denote  the  (2L  +l-rf)X(2L  +1)  matrix  constructed  (offline) 
such  that  4>(i )  defined  by 

*(*)  -  CflO  (5) 

satisfies  the  end  point  conditions4 


‘That  is,  $4i=l  if  k—  1,  while  $4 j=0  if  ky*l. 

*  It  is  easily  shown  that  the  matrix  C  has  full  rank  and  that  each  of  its  (2 L  +l-ff )  row  vectors  is 
determined  by  the  solution  to  a  Vandermonde  type  matrix  equation  [5,6[. 
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4>(* \o)  =  ^\T)  =  0,  i=0,l  .  .(n-1).  (®) 

This  construction  results  in  a  (2L+l-n)  dimensional  vector  valued  modulating  function  of 
order  n  in  which5  the  derivatives  p‘4>(t  )=<t )  have  the  representation  (cf.  (4)) 

(-l)V*(f)  =  CD^t),  i=0, 1,2..  (7) 

where  D  is  a  block  diagonal  matrix  defined  by 

D  -  w0  Dbg  [  0,  (  °j  J),(°2  g)  ..  (“i  ]  (8) 

and  D°  is  defined  as  the  identity  matrix. 

In  accordance  with  the  modulating  function  approach,  (3)  is  multiplied  by  <F(f )  and 
integrated  by  parts  n  times  while  noting  (6)  and  (7),  thereby  obtaining  the  vector  equation 

o  =  CY  +  £  £  CD*-"6"  z(0,k)ai(k)+  £  £  CD^ZUmU)  (9) 

+  EEE  cDi-'z(j,k)li(j,k) 

1=1;  =lt  =1 

where  ( Y  ,Z{j  ,k  )),  0</ <m  ,  0<*  </,  are  vectors  of  finite  Fourier  series  coefficients  of 
moments  of  the  data  defined  by 

T 

z(/,*)  -  /!“(<)]' do) 
0 

Y  =  Z( 0,1). 

These  vector  valued  functionals  of  moments  formed  from  the  data  can  be  calculated  using  an 
FFT  algorithm  as  discussed  in  [8j. 

After  introducing  some  notation,  a  rearrangment  of  (9)  leads  to  the  following  linear  equa¬ 
tion  in  the  coefficient  parameters: 

-CY  =  £  £  £n{m,l;i,j,k)${i,j,k)  (11) 

,=I;=U=1 

where  the  (2 L  +l-n  )X3  matrix  H{m  ,l;i,j,k)  and  3X1  vector  0(i ,j,k)  are  defined  by® 

ff(m  ,k)  =  CD  *_I  ±Df>'Z(0,k)lzU,0),ZU,k)  (12) 

m  I 

and 

=  fiiU) 
iiUt) 

respectively. 

4  Notice  that  4>(f )  is  also  a  modulating  function  of  order  n  for  any  n  <n  . 

*  The  reason  for  including  the  reciprocals  (l/m  ,1/f)  in  the  first  two  columns  of  (12)  is  to  account 
for  the  expansion  of  the  first  two  double  sums  in  (9)  to  a  triple  sum  so  that  (11)  and  (9)  are  equivalent. 
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Except  for  the  re-parametri2ation  and  notational  differences,  (11)  is  equivalent  to  the 
regression  equation  derived  in  [8]  and,  as  such,  can  be  used  as  the  starting  point  for  formulating 
a  least  squares  estimate  of  the  coefficient  parameters  under  the  proviso  that  the  structural 
parameters  (n,m,l)  are  known.  However,  the  emphasis  here  is  on  obtaining  values  (n0,m0,/0) 
which  are  correct  insofar  as  the  given  input-output  date  is  concerned.  For  this  purpose,  it  will 
be  assumed  that  (11)  contains  a  sufficient  number  of  algebraic  equations  to  permit  a  one-shot 
solution,  pending  some  nondegeneracy  conditions  on  the  data,  regardless  of  the  actual  values 
(n0,m0,/0)  as  long  as  (2)  is  satisfied.  Since  the  number  of  unknown  coefficient  parameters  is 
possibly  as  large  as  n(l  +m  +Tm)  and  since  the  number  of  rows  in  the  matrix  C  is  (2L+l-»), 
it  follows  that  L  must  be  chosen  in  the  general  nonzero  input  situation  to  satisfy  the  inequality: 

2L  +l-n  >  n(l  +  m  +Tm )  or  2  L  >  n  (TVm  +7fn+l  )-l.  (13) 

If  on  the  other  hand  the  input  is  zero  on  *0,7’),  i.e.,  u(t)=0,  0<t  <T ,  and  if  the  time 
interval  [0,T]  has  been  preceded  by  an  interval  during  which  all  the  system  modes  have  been 
excited,  then  (3)  reduces  to 

0  =  »(<)+  £  S«((*)p"1'f4"[»(OI* 

i=l4=l 

and  (11)  can  be  replaced  by  the  simpler  expression 

-CY  =  E  E H0(i,kMi,k )  (14) 

1=14=1 

where  the  (2 L  +l-n  )  dimensional  column  vector  HfJ[i,k )  and  the  scalar  90(i  ,k )  are  defined  by 

H0(i,k)  =  C0’~I+<“  Z(0,k)  and  90(i,k)  =  «,(*).  (15) 

In  this  special  situation  there  are  at  most  nT  unknown  coefficient  parameters  that  can  be  deter¬ 
mined  given  the  data  |0,y(<)l  on  [0, T’ ]  so  that  the  condition  on  L  can  be  replaced  by  the  less 
stringent  inequality 

2 L  >  W(T+l)-l.  (16) 

In  summary  for  this  section,  the  problem  of  structure  determination  for  the  class  of  poly¬ 
nomial  differential  systems  modeled  by  (3)  can  be  stated  as  finding  the  smallest  triple  indexed 
“discrete  time”  integers  (n0,m0lI0)  and  “control  vector”  triple  indexed  sequence  9{i,j,k)  such 
that  (11)  is  satisfied  given  the  discrete  indexed  “system  matrix”  ,/ ,k)  and  the  output 

vector  CY .  This  problem  will  have  a  solution  only  if  L  has  been  chosen  to  satisfy  (13),  or  less 
restrictively  (18)  if  the  zero  input  prevails  on  [0,7’J,  and  only  if  the  vector  CY  is  reachable  for 
some  (n,m,l)  satisfying  (2),  a  condition  which  reflects  back  on  the  input-output  data 
[u  (I  ).lf(* )]  on  [0,  T ]  as  will  be  seen  below. 

3.  MINIMUM  NORM-TIME  OPTIMAL  SOLUTION 

As  formulated  above,  structure  determination  is  a  simplified  version  of  a  discrete  system 
optimal  control  problem  considered  by  Sarachik  and  Kranc  [9],  although  the  discrete  time  is  tri¬ 
ple  indexed.  Following  the  analysis  used  in  [9J,  consider  the  linear  functional  derived  from  (11): 

-X'CK  =  E  E  EV//(m,/;.' ,/,*W ,/,*)  (17) 

i  =iy =1*  =1 

where  X  is  an  arbitrary  (2 L  +1-R  )  dimensional  column  vector.  The  right  hand  side  of  (17)  will 
be  viewed  as  an  inner  product  for  fixed  values  of  (n  ,m ,/)  by  defining  the  3X1  vector  function 
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K{i,j,k)*s 

K(iJ,k)  =  //'(m  )X. 

Using  (18)  and  the  inner  product  notation  <  ,  >,  (17)  is  represented  by 

-\'CY  =  <K,0>. 


Taking  absolute  values  of  both  sides  of  (19)  and  applying  the  Schwarz  inequality  to  the  right 
hand  side: 


WCY\  <  I  A  im 


where  the  norm  on  $  is  defined  by 


‘  •  m  l  l* 

M  =  EE  E *'(•*,/.* W,/,*) 

•=i>  =u  =i 


and  the  norm  on  K  is  similarly  defined. 

From  the  inequality  (20)  it  can  be  deduced  that  if  the  problem  has  a  solution  for  an  arbi¬ 
trary  vector  CY ,  it  must  be  the  case  that 

M  >  U'l1 


where  K* (i  ,j  ,k)=H'(  m,/;i  ,j,k))f  and  X*  is  that  value  of  X  which  solves  the  minimization 


problem 


Min  I  K  |  =  X'E  E  E  H'(tn,l;»,j,k)H{m,l,iJ,k)\ 

i=iy=u=i 


subject  to 


-x'or  =  i. 


The  latter  problem  is  equivalent  to  minimizing  the  quadratic  form 

I  K[2  =  X'  W(n  ,m  ,1  )X  +  n(l+\' CY)  (22) 

where  /i  is  a  scalar  Lagrange  multiplier  and  W(n  ,m  ,/)  is  the  Gram  matrix  defined  by7 


W(n,m,l)  =  £  E  E  H{m,l  ,i,j,k  )H'{m,l  ;i,j,k ). 

i=l;=l*=l 


Assuming  the  inverse  of  lF(n,m,/)  exists  for  some  (n,m,f),  the  minimizing  value  of  X  subject 
to  the  indicated  constraint  is  found  to  be 


X*  =  W~l{n  ,m  ,l)CY 
~  Y'C’W-l{n,m,l)CY 


If  a  solution  to  the  above  minimization  problem  exists  for  some  (n0,m0,fo)  satisfying  (2),  a 
solution  exists  for  all  (n  ,m  ,/)  satifying  the  following 

n0<n  <n,  m0<m  <m ,  l0<l  <T.  (24) 

This  is  due  to  the  monotonicity  property  of  the  Gram  matrix  (23)  in  that  using  matrix  inequali¬ 
ties  in  the  sense  of  symmetric  positive  semidefiniteness: 


7  Note  that  W(n  ,m,l)  can  be  interpreted  as  an  output  controllability  Grammian  for  the  problem 
at  hand. 


W(n0,m0,l0)  <  W(n  ,m  ,/ ) 

holds  for  all  (n  ,m  ,/)  satisfying  (24).  This  coupled  with  the  fact  that  the  norm  (21)  is  a  mono- 
tonically  increasing  function  of  ( n,m,l )  allows  for  the  conclusion  that  structure  determination 
for  the  model  (3)  given  the  data  at  hand  corresponds  to  the  smallest  triple  of  integers  which 
satisfy  the  determinant  inequality: 

det  tV(n0-l,ro0-l,/0-l)  <  e0  (25) 

together  with  one  or  more  of  the  following  inequalities: 

det  W(n0,m0,i0)  >  (0  or  det  W'(n0-l,m0,/o)  >  <o 

det  W(n0-l,ro0-l,/0)  >  «o  or  det  W^no-l.mo./o-l)  >  fo 
det  W(n0lm0-l,l0-l)  >  t0  or  det  H^Oo, mo-1, f0)  >  fo 
det  ir(n0,m0,/0-l)  >  e0 

where  f0  is  a  positive  parameter  selected  by  the  user  to  discriminate  between  the  theoretical  zero 
value  and  small  threshold  effects  of  noise,  roundoff  errors,  and  the  like. 

Nothing  is  claimed  about  the  existence  or  uniqueness  of  such  a  triple  (n0,m0,/0)  satisfying 
(25).  About  all  that  can  be  claimed  at  this  point  is  that  if  a  minimizing  triple  can  be  found 
satisfying  (25),  then  that  triple  (n0,m0,/0)  solves  the  structure  determination  problem  relative  to 
the  data  at  hand.  More  analysis  is  needed  to  find  conditions  on  the  input  u  (t )  which  would 
guarantee  the  existence  of  a  solution.  The  problem  is  complicated  by  the  fact  that  the  modulat¬ 
ing  functions  act  as  a  filter  on  the  data  and  that  the  various  polynomials  arising  from  the 
differential  operators  in  (3)  may  have  some  common  factors,  i.e.,  they  may  not  be  coprime. 
Nonetheless,  (25)  is  a  relatively  easy  test  to  perform  in  consideration  of  the  complexity  of  struc¬ 
ture  determination  problems,  and  any  minimizing  solution  to  (25)  is  a  likely  candidate  for  the 
structure  parameters  pending  additional  tests  involving  a  variety  of  inputs,  initial  conditions, 
etc. 

Assuming  the  existence  and  determination  of  a  minimizing  triple  (n0,m0,/0)  satisfying  (25), 
the  system  coefficient  parameters  for  (3)  can  be  obtained  using  the  least  squares  procedure  of  [8], 
or  if  sufficient  confidence  exists  for  the  one-shot  data  at  hand,  the  conditions  for  equality  in  the 
Schwarz  inequality  (20)  yield  the  parameters  directly  as 

**(»,/.*)  =  ,j  ,k)W-\n0,m0,l0)CY  (26) 

l<i<n0,  1  < ;  < m 0,  1  <*  < /0- 

In  the  special  case  of  the  zero  input  on  [0,  Z* ],  W(n  ,m ,/)  is  to  be  replaced  by 

wv«.0  =  £  £//(>(»•, A  MM)  (27) 

•  =1  i  =1 

where  H0(i,k)  is  defined  by  (15).  Assuming  existence  and  determination  of  a  minimizing 
integer  pair  (n0,/0)  satisfying 

det  JVo(n<rl,/0-l)  <  (0  (28) 

together  with  one  or  more  of  the  following  inequalities: 


(29) 


j 

»' 

i 

i 

» 


i 
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det  Wo(n0,/0)  >  <o  or  det  W^o^o-l)  >  <o 

det  W0(n0~l,l0)  >  (0, 
the  corresponding  relation  to  (28)  in  this  special  situation  is 

*£(•'.*)  =  **(*)  =  -H()(i,k)W0-i(no,l0)CY 

l<i<n0, 


4.  CHOICE  OF  (L,T)  AND  NOISE  CONSIDERATIONS 

Suppose  the  system  bandwidth  wj  is  approximately  known.  Although  the  value  of  wj  is 
likely  to  be  input  dependent  for  a  nonlinear  system,  it  seems  reasonable  to  presume  that  such  a 
value  may  be  available  from  practical  considerations.  Clearly  the  highest  frequency  Lu0  in  the 
modulating  functions  should  be  chosen  comparable  to  the  system  bandwidth  else  higher  fre¬ 
quence  noise  will  be  present  in  the  vectors  Z(j  ,k),  i.e. 

Lu0  =  L  2ff/  T  ~  w*  .  (30) 

Assuming  the  inequality  (13)  is  approximately  satisfied  with  equality,  it  then  follows  that  a 
guideline  for  choosing  T  is 

T  *=»  (n(/+m+7rn+l)-l)jr/u;j  .  (31) 

In  the  special  case  of  a  zero  input  on  [0,T],  the  analogous  relation  for  a  guideline  in  choos¬ 
ing  T  =  T0  is 

r0~  (n(r+l)-l)ff/wj .  (32) 


5.  CONCLUSIONS 

The  determinant  inequalities  derived  in  (25)  and  (28)  indicate  ways  of  testing  for  the  struc¬ 
tural  parameters  of  the  model  (3),  given  upper  bounds  for  the  indexes  and  assuming  the  system 
belongs  to  this  class  of  nonlinear  models.  Although  relatively  simple  in  appearance,  this  simpli¬ 
city  undoubtedly  belies  a  variety  of  the  difficulties  that  might  be  encountered  in  a  practical 
situation.  Not  only  will  the  presence  of  significant  amounts  of  noise  invalidate  the  conclusions, 
but  there  are  potential  difficulties  resulting  from  the  projections  on  the  data  due  to  the  modu¬ 
lating  functions  and  the  possibility  of  common  factors  in  the  polynomials  corresponding  to  the 
differential  operators  in  (3).  The  noise  problem  can  be  mitigated  to  some  extent  by  a  careful 
choice  in  modulating  function  frequencies  as  explained  earlier.  In  fact,  there  is  no  need  to  pick 
the  first  L  commensurable  sinusoids;  rather,  the  basis  set  (sinfw0I  ,cos/w0t },  0<t  <T , 
w0=2jt/  T ,  l£(l „  .  .  lL)  will  suffice,  where  (/,,  .  .  lL)  is  any  suitably  chosen  integer  set.  The 
possibility  of  modes  in  the  data  which  are  orthogonal  to  one  or  more  of  the  modulating  func¬ 
tions  is  practically  unimportant  if  the  number  of  such  modulating  functions  is  sufficiently  large, 
i.e.,  L  is  chosen  sufficiently  large,  to  encompass  all  the  essential  information  in  the  data.  How¬ 
ever,  such  a  choice  may  then  require  a  very  long  time  interval  [0,r]  or  high  frequency  noise 
problems  may  be  encountered,  etc.  Further  analysis  will  be  needed  to  resolve  these  questions. 
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